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Q^. Abstract 

w\| In this work, we are investigating the problem of integrability of Bianchi class A 

cosmological models. This class of systems is reduced to the form of Hamiltonian 
systems with exponential potential forms. 
r> I The dynamics of Bianchi class A models is investigated through the Euler- 

C^ ■ Lagrange equations and geodesic equations in the Jacobi metric. On this basis, 

we have come to some general conclusions concerning the evolution of the volume 
function of 3-space of constant time. The formal and general form of this function 
has been found. It can serve as a controller during numerical calculations of the 
dynamics of cosmological models. 

The integrability of cosmological models is also discussed from the points of 
view of different integrability criterions. We show that dimension of phase space of 
Bianchi class A Hamiltonian systems can be reduced by two. We prove vector field 
of the reduced system is polynomial and it does not admit any analytic, or even 
formal first integral. 



1 Introduction 

We shall investigate the dynamics of the most interesting group of homogeneous Bianchi 
class A cosmological models which is described by the natural Lagrangean function 

^ = \9.prq^-V{q)=T-V{q) = 

1 ^ ding' din q^ ^ (o \- i J ^ 2 2\ f,. 

where qi ~ A\{i = 1, 2, 3) are three squared scale factors Ai for diagonal class A Bianchi 
models; different Bianchi types correspond to different choices of rij G {—1,0,1}, i = 
1, 2, 3; a dot denotes differentiation with respect to cosmological time t. Logarithmic time 
T is related with cosmological time t by 

, dt dt 

dr 



^qiq2q3y/6 VoliVfS ' 

Bogoyavlensky |l|, proves an important property of system (|l|), namely the existence 
of the monotonic function F with the following form 

F = Jt(^'^"^')'^' = Jt^^^^'^ (2) 

such that 

dt ~ 
Function (^ is invariant with respect to the scaling transformations and it has the sense of 
the speed of change of the average radius of the universe. Function \F\ along any solution 
decreases from infinity to zero in such a way that F = is reached at the moment of 
maximal expansion, and |F| = cx) corresponds to the initial singularity. The existence of 
function F allows us to define what we call the early stage of the evolution of the universe 
as 

F>1 . 

The importance of this function for numerical integration of B(IX) models has been 
pointed out in the work 0. The authors have used the Rauchaudhuri equations to show 
the property of upper- convexity of function (VolM^)(t) which means that this function 
does not possess a local minimum (where F = and F > 0) and may possess not more 



than one maximum. If F < 0, the volume function (VolM'^)(t) shrinks; whereas if F > 
it expands. Both of the processes take place in the same region of the phase space {p, q) 
but with reverse directions of time. In the phase space (p, q) function F has the following 
form 

{q\Wf" ^ . dL 

P = 5 Z^P*^' P^=d^^ 



i=l 

and 

lU ~ 9 



T.P^<f] -6V^ 



^4 = 1 



where pi are momenta conjugated with generalized coordinates g*. 

In work [Q , function F was used for controlling the quality of numerical integrations of 
the B(IX) model. In this model, the scale factors oscillate in a neighborhood of the initial 
and final singularity. The function (VolM^)(t), obviously, does not possess the analogous 
property g, |]. 

In the present work we give some general and formal expressions for the function 
(VolM^)(t). It can be used as a tool for studying the B(IX) models. Let us note recent 
important results of Cushman and Sniatycki concerning chaos in the B(IX) system. 
They proved that existence of a monotonic function F excludes the possibility of recur- 
rence in the system and, thus, any form of standard deterministic chaos in the system. 
This illuminates previous negative results and shows that for a study of this system we 
have to use non-conventional methods. 

Several authors tested if the last model passes the standard Painleve integrability 
test (in the form of the ARS algorithm |1^). First results of Contopoulos [Q shows that 
B(IX) model passes this test. Next, this paper was revised P, however, without any strict 
conclusions concerning integrability. It was stated also that this model passes Ziglin's test 
(see p3|, |2^). More careful Painleve analysis was done by Latifi et al. in |^. They show 
that B(IX) model does not passes the so called perturbative Painleve test. Authors of 
this paper suggest the existence of 'some chaotic regimes' in the system. 

The above remarks show that the notion 'chaos' has unclear status when dynamical 
systems arising from the general relativity and cosmology are studied. Moreover, for B(IX) 
model discrete dynamics defined in |]^ shows strong ergodic properties, however, this 
'chaotic behaviour' seems to be absent (or hidden) in the continuous dynamics. Moreover, 
the standard criteria of detection of chaos (Lyapunov characteristic exponents — LCE) are 
not invariant with respect to the time reparametrization and transformation of phase 
variables whereas existence of first integrals is an invariant property of the system. It is 
also important to note that the non-zero LCE can be used as an indicator of chaos only 



when the motion take place in a compact invariant subset of the phase space, but it is not 
true for the B(IX) dynamical system. All these facts motivated us to study the problem 
of integrability of the investigated models. 

The non-integr ability of the system is a weaker property than chaos (in the sense of 
the deterministic chaos) but better described and understood. The authors believe that 
investigation of non- integrability in B(IX) models can contribute to a better understanding 
of chaos in cosmological models. Here we show that the Bianchi class A Hamiltonian 
system are not completely integrable in the sense of Birkhoff. This conclusion is weak as 
the negative answer to the question about algebraic complete integrability of B(IX) (see 
[pl[]). In order to obtain stronger result we reduce the dimension of the phase space by 



two. We show that the reduced system is polynomial and, what is most important, it 
does not admit any analytical, or even formal, first integral. 

In cosmological models chaos, if properly defined and present, has some hidden char- 
acter. The basic indicator of chaos in these models, the LCE, depends on the choice of 
the time parametrization. In the logarithmic time r, nearby trajectories diverge linearly 
whereas in other time parametrizations they will diverge exponentially which is charac- 
teristic for chaotic systems. The fact that the rates of separation of nearby trajectories 
depend on the clock used is obvious. The problem is in invariant choices of the time 
parameter for the invariant chaos detection. Such a role is played by Maupertuis clock 
(time parameter s is such that ^ = 2\E — V\, where E is the total energy of the system, 
V is its potential and r is mechanical time). 

Our point of view is such that the LCE, when used in general relativity, should be 
defined in an invariant way. Then the results could be interpreted in a different time 
parametrizations. The Bianchi IX model is 'chaotic' in the parameter s (LCE is positive), 
but, after transition to the parameter r, nearby trajectories diverge linearly in such a 
way as integrable systems. This phenomenon is called the hidden chaos. Let us note that 
the existence of the first integral of an autonomous system is an invariant property (with 
respect to time reparametrization and to transformation of phase variables). 

In general relativity and cosmology, the problem of non-integrability or chaos is not 
only very subtle but also is strictly connected with the invariant description. One must 
be very careful detecting integrability in B(IX) dynamics. The problem whether chaos in 
the gauge theory is a physical phenomenon is, generally, open. 



2 The dynamics of Bianchi class A models from the 
Euler-Lagrange equations 

The Hainiltonian function for the system (1) has the following form 

^=^rtaP/3 + n9), (3) 



where 

/ — (g^)^ q^q"^ q^q^ 



Pa = gafsq^, 



g"^ = 2 



q^q^ —{q^Y q^q^ 

\ q^q^ q^q^ —{q^Y 



\ i<j i=i 

V{q) is the potential function. The obtained Hamiltonian system is considered only on 
an invariant set of the phase space defined by the zero level of the Hamiltonian @, i.e., 

n = 0. (4) 

The Euler-Lagrange equations in time r have the following form 

dr' ^^^^dr dT~^ dqP' ^^' 

where Christoffel symbols F^ are connected with the metric defined by the kinetic energy 

T = Igafs^q^ = Ig'-^PaPp. 

Equations (^) after transformations to a new time parameter s, called Maupertuis time, 
take the form of geodesic equations for the Jacobi metric 

^„^ = 2\E - V\gai3 = 2Wgaf3, 

that is, 

d\- ^„ dq^ dq^ 
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where a hat denotes that respective quantities are calculated with respect to the Jacobi 
metric. Christoffel symbols calculated from g and g metrics are connected by relations 

where 

$ = - In 2W. 
2 

Let us note that the kinetic energy form does not depend on the Bianchi type models 
characterized by the set {ni,n2,n3}. The only non-vanishing Christoffel symbols are 

^11 = "qi' ■'"22 = --^> Tgg = - — . (8) 

After substitution (H), system (|) takes the form 

^^ - (^1^^^) = K)'(^'T + (n.riq'Y - {n,f{q^f - 2n,n.g^g^ (9) 

where {2, j, k} G 6*3, and S^ denotes the set of even permutations of {1, 2, 3}. 
The change of variables 

q' = e^\ 2 = 1,2,3, (10) 



transforms the above equations to the following form 

^n^fe^Q^ + i^nufe'^" - {n^fe^^' - 2n,n,e«^+«\ (11) 



cPQ' _ 



where {i, j, k} G S'3. 

The system ([TT|) is satisfied on the Hamiltonian constraint 7i = which is equivalent 
to the condition of normalization of the tangent vector to the trajectory m* = dq^ /ds, that 
is, 

\\uf = '^^9ai3^^ = -sgnl^, (12) 

or 

dq" dq^ 
9af3-^-^ = 2l^sgnK 
dr dr 



In terms of variables Q*, the constraint condition is equivalent to 

dr dr 



J^ dO' dO^ , , 

E^^ = -8^- (13) 



Adding the sides of equations ([TI|) , we obtain the following formula 

dr2 



Equations (pJ]) , after time reparametrization r ^* s = s{t)^ take the form of geodesic 
equations. From (^) we obtain 

2n» ^n* f)\M Ani 



^^'^ + ^1^^ = ^^^■^'"'"'^ + ^^^)'''''' " ^^^)'''''' " 2n,n,e«^+«'= 



where {i, j, h) G S'3. 

The problem of investigation of Lagrange systems with the indefinite kinetic energy 
form is open. First steps to investigate such systems have been done in [§]. In the 
terminology used in ^ our system is a special case, the so called non-classical simple 
mechanical system. As it was established, these systems have the following fundamental 
property. A trajectory of the system can pass through the set dD = {g : i? — V^ = 0}. 
During this passage the vector tangent to trajectory changes the cone sector defined by 
the kinetic energy form: gaf^iQo)^"^'^ = where ^" = dq°'/ds, go ^ dD. In our case the 
signature of ga/3 is Lorentzian, i.e., (— , +,...,+) (for details see 0). 

In generic situations {rii 7^ for i = 1, 2, 3) which include BVIII and BIX models 
(Mixmaster models), there are analytical and numerical arguments that the function of 
sign of the potential for a typical trajectory is an infinite subsequence that is a one sided 
cut of the following double infinite sequence (see |jlO[) 

sgnV^ = {...,+1,0,-1,0,-1,...}. (15) 

If we assume that the subsequence (|TB|) is finite, then our system reaches the state V = 
{W = in general case) which corresponds to Kasner solutions, finite number of times. 
During the Kasner epoch the information about the localization of the point on the interval 



of normal separation (modulo initial localization) grows e-times [TT| . If the subsequence 
(|15D is finite it means that after n epochs n bytes (i.e., finite number of bytes) of infor- 
mation have been lost whereas we know that our system is chaotic (the loss of infinite 
information is required for chaos). Let us notice that when the system goes asymptotically 
to the boundary sets W = then it is asymptotically free. 



3 Properties of the function of volume of the con- 
stant time 3-space 

Equation ([T^) implies that 

1 d^ 

-— ln(VolM3) = ~2V = 2T. (16) 

The above relation means that, in a generic case (Vinj ^ 0), there is an infinite number 
of intervals in which function In(VolM^) is subsequently convex up and down. These 
intervals are separated by an infinite number of infiexion points (which corresponds to 
^ = 0) in the diagrams of the function ln(VolM^)(r) and lie on the lines ln(VolM^)(r) = 
±r + C. The additional information we have about the B(IX) model is that this model 
has the initial and final singularities. In the following paragraphs of the work we shall 
concentrate on the B(IX) models. 

From the (0, 0) components of the Einstein equation for the BIX case, we obtain 
that function ln(VolM^)(r) cannot possess a local minimum but it can possess a single 
maximum. The above property suggests that for a typical trajectory in this model the 
qualitative diagrams of the function ln(VolM'^)(r) look like in Figure 1. 

Fig. 1. 



After integrating over r the both sides of (JT^ and assuming that in the moment of 
maximal expansion r = tq, we obtain 

ln(VolM3)(r) oc e^^^J^o «Wsgn(-v).*^ ^^^^ 

where we choose Ci = 1 in the phase of expansion and Ci = —1 in the phase in contraction 
of the volume function (if V^ = 0, ln(VolM^)(r) oc e^"^). Finally, for any model which 
describes the evolution of the volume function 



ln(VolM=^)(r) oc e"(^+<^>) ^r^^ e"e<'> oc e\ {U 

{s) = - r s{t)sgn{-V)dt, 

r Jrn 



'TO 

where we assume that the average value of s(r) on the interval (r, tq) exists as r 



and it is finite. From the formula ( |18D it immediately yields that in a neighborhood of the 
initial singularity (r —>■ — oo) the volume function changes exactly as in Kasner's models. 
The second observation is as follows: the volume function does not oscillate around the 
equilibrium positions ln(VolM^)(r) = but oscillates around Kasner's solution. In other 
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words, Kasner's solution plays the role analogical to the equilibrium positions in the small 
oscillation approximation. 

From ( ]T7| ) one can obtain the following relations between a natural parameter s defined 
along geodesies (Maupertuis time) and the volume function ln(VolM^)(r) 

s = 4ln(VolM^)(r) for V<0, 

at 

s =-— ln(VolM3)(r) for V > 0. 

at 

The zero value of the parameter s corresponds to the moment r = tq. The relations 
between the parameter s, function F and the scalar expansion function 

e = -^in(VolM=^)(t), 

are as follows 

s{t) = ±6[(VolM=^)(r)]^/=^F(r), 
s{t) = ±[(VolM3)(r)]e(t(r)), 

where plus and minus sign correspond to V < 0, and V^ > 0, respectively. 

From the above, we can conclude that "near the singularity" is equivalent to s ^ 
(Vo1M3)V3^ i.e., s > 0. 

Formula (|TBp implies that the characteristic time after which (VolM'^)(r) grows e- 
times, i.e., (VolM^)(r) oc e'^/'^'^'"'"' has the following form 

Tchar = (1 + {s))-\ (19) 

This characteristic time is finite if the average (s) exists. 

For a typical trajectory the function s{t) is shown in Fig. 2. 

4 The Bianchi class A models as systems with expo- 
nential potentials and Its Algebraic Non-integrability 

After introducing the new variables Q* and using the definition (10), the Lagrange system 
(1) can be transformed to a Hamiltonian one. The Hamilton function for this system 



takes the following form 

i<j i=j \ i<j i=j 

= ^^"Vp/3 + V^(Q"), (20) 

where pi = \{Q^ + Q'') for {i,j, k} G 5*3. The Hamilton function pO| ) is a special case of 



the Hamiltonian for the so called perturbed Toda lattice [11 



5 Analysis of the integrability of B(IX) model 

There are several definitions of integrability. Generally, integrability means that the 
system under consideration possesses a large enough number of first integrals. Necessarily 
we have to specify the class of functions that contains these first integrals as well as 
to define the domain of their definition. Let us note here that there are examples of 
Hamiltonian systems that possess an first integral of class C" but does not possess an 
integral of class C^ with (5 > a, ioi a, P = 1,2, . . .oo,uj (see |I8|). 

It is well known also that every system of ra differential autonomous equations is locally 
integrable — in a neighborhood of every nonsingular point (where the right hand sides do 
not vanish) it possesses n — 1 first integrals. Thus, nontrivial problems are non-local or 
concern the existence of integrals in a neighborhood of equilibria points. 

It is very difficult to prove (non) integrability of a given set of differential equations. 
One way to simplify the problem is to restrict a class of function where we look for 
integrals. As an illustration of this approach, let us consider the Birkhoff integrability 
(see [|l^) of the Bianchi class A system in the form (|20|) . This system belongs to the wide 
class of Hamiltonian systems in R^" equipped with the standard symplectic structure and 
are given by the following Hamiltonian function 

T^ = ljiP,p)+ J2 Vm.ew{cm,q), (21) 

meM 

where 



n 



{P,P) = J2 (^'^PiPj^ i'^m, q) = Y1 '^"i,q\ m = (mi, . . . , m„) G Z", 
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(a*-'), Vm and Cm are constant; A1 is a finite subset of Z" 

7W = {m e Z" I t;^ ^ 0}. 
We look for integrals tliat are polynomials with respect to p, i.e., 

fiq,p) = E /fcp^ 

where 

n 

Xi = {meZl\\m\<l}, \m\=Y,mu p"" = pT' ■ ■ -PT^ m E Zl, 

i=l 

and coefficients fk have the form of infinite series of exponents 

mGZ" 

Here Z+ denotes non-negative integers. We say that the system (^) is Birkhoff integrable 
if it possesses n independent integrals of the prescribed form (see note of Ziglin [^ about 
modification of the original definition of Kozlov). We order elements of A4 with respect 
to the lexicogrphic order and denote by a its maximal element and by (3 the maximal 
element of M that is not colinearly with a. Then, according to Theorem 3 from |^ if 

k{a, a) + {a, (3) ^ 0, for all k e Z+ (22) 

then the Hamiltonian system ( p^ is not integrable in the sense of Birkhoff. We immedi- 
ately have 

Theorem 1 A generic case of the Bianchi class A system given by Hamiltonian function 
^l\ ) with Tii j^ for i = 1,2,3 is not integrable in the sense of Birkhoff. 

Proof. For the Hamiltonian ([2^) we have 

A^ = {(1, 1, 0), (1, 0, 1), (1, 0, 0), (0, 1, 1), (0, 1, 0), (0, 0, 1)}, 

and thus a = (1, 1, 0) and /3 = (1, 0, 1). Metric (a*-') has the form 

" -1 1 1 

{a'^)=2 1 -1 1 

_ 1 1-1 

and thus we have 

k{a,a) + {a, P) = A ^ 0, 

and this finishes the proof. □ 

Let us remark that in a case when one of n, is equal zero then k{a, a) + {a, (3) = for 
all k. In such a case the system has one additional integral, namely pi. 

11 



6 Reduction and Non-integrability of B(IX) model 

The results obtained in the previous section are weak. There are two reason of it. First, 
we asked about the complete integrabihty of the system. However, the system under 
investigation can have only one additional integral. The most important is the fact the 
we pose our question for the system defined on R^ although we are interesting only in 
the system on five dimensional manifold defined by the level 7Y = 0. One can imagine a 
system that is not globally integrable although it is integrable on a one prescribed energy 
surface. 

In this section we want to study the B(IX) Hamiltonian system just only the level 

n = o. 

During investigation of a dynamical system we usually try to lower its dimension 
making use of its first integrals and symmetries. For Hamiltonian system (^) we know 
only one first integral — Hamiltonian. Thus, using it we can potentially can make reduce 
the dimension of the system by one however we loose the polynomial form of the system. 
However we do not want this side effect of reduction because it excludes possibilities of 
applications algebraic tools for study non-integrabihty. 

In this section we show how to reduce the dimension of the phase space by two and 
preserving the polynomial form of the considered vector field. In what follows consider 
her the case of B(IX) {rii = n2 = n^ = 1). 

First, we transform the Hamiltonian vector field corresponding to Hamiltonian (^), to 
a homogeneous polynomial form of degree two. To this end let us put 

yi = q\ Zi = — « = 1,2,3, (23) 

then equation of motion will have the form 

iji = yiZi, Zi = {yj - y^f - yf, {z, j, k} e S^. (24) 

This system has the first integral corresponding to the Hamiltonian (Q). It is has the form 
H = Z1Z2 + ZiZs + Z2ZS -yl + 2yiy2 -yl + 2yiy^ + 2?/2l/3 - yl- (25) 

We make change of variables 

'Wi = yi + 1/2, W2 = yi- 2/2, w^3 = ?/3, (26) 
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and we leave Zi unchanged. In new variables system (|2^) has the form 

wi = lzi{wi + W2) + \z2{wi -W2), 

W2 = lzi{wi + W2) - lz2{Wi - W2), 

W3 = Z3W3, 

Zl = {W3 -Wi){w2 + W3), 
Z2 = {W3 -W2){W3 -Wi), 
Z3 = {W3 + W2){W2 -tys), 

and the first integral (^) is transformed to the following form 

H = Z1Z2 + Z1Z3 + Z2Z3 - wl + 2wiW3 - W3. 
Now, we introduce new variables 

Zl Z2 Z3 W2 Wi 
Ml = , U2 = , U3 = , M4 = , U5 = , Ue = W3. 

W3 W3 W3 W3 W3 

After this transformation we obtain the following system 

ill = Ug[{1 + M4)(l - Mg) - M1M3], 
U2 = M6[(l - M4)(l - M5) - U2U3], 

U3 = ueiul -uj-l) 

ilA = IU6[U4,{UI +U2- 2U3) + M5(mi - U2)] 

U5 = lue[u4,{ui - U2) + u^{ui +u2- 2u3)u5] 



(27) 



(2^ 



(29) 



U6 = U3UQ, 



(30) 



with the first integral 



H = ul{uiU2 + U1U3 + U2U3 - ul + 2m5 - 1) 



(31) 
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Now, we make use the fact that B(IX) model is considered only on the level H = 0. From 
equation if = we find u^ as a function of (mi, U2, M3, M4): 



M5 



1 



1 + U4 - U1U2 - U1U3 - U2U3) 



thus we can eliminate this variable from the right hand sides of (pOl). Moreover if we 
change the independent variable according to the rule 



d 
di 



uq d 
Yds'' 



(note that Uq > 0) than the first four equations in ( pOf ) separate from the last two. Thus, 
we finally obtained the following close system describing the dynamic of the B(IX) model: 



Ui = {1+ M4)[l + U1U2 + U^iui + U2) - UJ] - 2UiM3, 



U2 



1 - M4)[l + U1U2 + U^iui + U2) - UJ] - 2u2Us, 



U 



■4J- 



(32) 



M3 = 2{ul -ul~l) 

U4 = U4{Ui +U2- 2Us) + |(mi - U2)[l - U1U2 - Us{ui + U2) 

This system will be called the reduced B(IX) system. We consider this system in C^ 
Theorem 2 The reduced B(IX) system does not have a non-trivial analytic first integral. 

Our theorem will be a consequence of the following lemma. 

Lemma 1 Consider a system of differential equation 

X = fix), /(O) = 0, fix) = (/i(x), Idots, fnix)) X e C", (33) 

with analytic right hand side, with 

fix) = Ax + Oi\x\'^), (34) 

where matrix A has eigenvalues Aj G C, z = 1, . . . , ra. If the system possesses an analytical 
first integral F then there exist non-negative integers ii, . . . ,in such that 



fc=i 



0, 



Y.ik>0 
fc=i 



(35) 
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Let us assume that an analytic fist integral exist and that condition ( p5D is not satisfied. 
We represent the first integral in the following form 

oo 

F = J2Fi, Ffc^O, k>l, 

l=k 



where Fi is a homogenous form of degree / 






' -^-K 



ik G Z+, fc = 1, Idots, n. 



From the equation 



5:/,(x)9,F=o, 



we conclude that the form Fk is a first integral of the system x = Ax, i.e., 



(36) 



Y,l,ix)d,Fk = 0, kix) = Y.A. 



ijXj 



(37) 



If matrix A is diagonalizable then we can assume that li{x) = XiXi, and then equation 



(|37D reads 



E 



<k) 



iiA \-in=k 

This equation implies that 



J2^i^i 



.1=1 



x\ 



•x!," =0. 



(38) 



<k) 



ll,...,lr, 



J2^i^i 



.1=1 



for all {ii,...,in) eZl,^ii = k, 



1=1 



<k) 



Because, F^ ^ there exit indices (^l, . . . , ^„) such that F^^^ _ ^^ 7^ and that for such 
indices we have 

n 

1=1 

Contradiction with our assumption prove the Lemma for the case of a diagonalizable 
matrix A. In the case of non-diagonalizable matrix the prove is only technically more 
difficult. See |^ and especially ||2^ where this approach is generalized. 

Let us remark the the above Lemma is also true if we assume the first integral is a 
formal power series. 
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To prove our theorem let us notice that for the reduced B(IX) model point z = 
{—i, —i, i, 0) is an equilibrium point and the matrix of the linearized system is diagonaliz- 
able and has possesses the eigenvalues {—2i, —2i, —ii, —ii). For this eigenvalues condition 
( |35[) cannot be satisfied and imply that the system does not have an analytic first inte- 
gral. In fact we prove more, namely, the system does not have a first integral that can be 
expanded around the point z into a formal power series. 

7 Conclusions 

The particular integrable subclasses of the Bianchi models play an important role in the 
analysis of the dynamics of the cosmological models. To illustrate this fact, let us con- 
sider the phase space of the solutions of Bianchi models in the Bogoyavlensky approach 
[0]. In the Bogoyavlensky methods of investigations of the corresponding dynamical sys- 
tems, we glue to the phase-space the boundary A onto which the system prolongs almost 
everywhere. The systems on the boundary A can be integrated and, in this way, we can 
study basic properties of the trajectories near the singularity. From the existence of the 
monotonic function F, we obtain that in the generic situation (Vi Ui ^ 0) the trajectories 
of the Bianchi class A models close up to the boundary F as F <^ — 1. 

Now, the trajectories move along the corresponding ones lying on the boundary. All 
the trajectories are the separatrices of the critical points. At last the trajectories reach 
the neighborhood of critical points K (corresponding to the Kasner asymptotics of the 
space-time metric) and they begin to move along their separatrices. The corresponding 
space-time metric for the mixmaster models is the BKL approximation [^. 

In this way the chaotic systems (so non-integrable) in Bogoyavlensky' s approach can 
be well approximated by an integrable system. This feature of so surprisingly good ap- 
proximation is not so far understood completely. 

In Bianchi models with chaos we have the infinite series of Kasner epochs (it would 
be good if we had a precise and exact proof of this fact) and these models do not exactly 
admit the Kasner asymptotics. 

It may be the case when we consider the higher-dimensional generalization of the 
mixmaster cosmological models (or models with a massless scalar field). These systems 
admit exact Kasner asymptotics, so, according to our theorem they will be integrated. 



In the work [jT6[ it has been proved that the BKL approximation is true not only for 
BVIII and BIX, but it exists for all Bianchi types (exceptionally BI and BV) with the 
movement of matter (n" ^ 6^). 

In the models I, V and also III, VII and the type II model, with the restrictions on 
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the velocity ui = U2 = u^ = 0, it is known that the Kasner solution is a general one near 



the initial singularity [|T6|. From the above, we can conclude that all these models are 
integrable. 

The facts that the BKL approximation represents a typical state of the metric in a very 
early state of the evolution has a very simple interpretation. Let for F = Fi -C — 1 there 
be some distribution of initial conditions (e.g. homogeneous). Approaching the initial 
singularity this distribution for F = F2, where F2 < Fi, transforms to the corresponding 
one concentrated in the neighborhood of the critical points. These critical points have the 
separatrices which move towards the physical region of the phase space. During the motion 
along such a separatrix the space-time metric is described by the BKL approximation. 

From the fact that F — >■ — 00 near the singularity, we can conclude the existence of 
a fundamental property of the system - the property of concentration of the trajectories 
near the boundary A. 

To finish with, we would like to make a certain suggestion more philosophical or 
methodological in character. Among the Bianchi class A models, the most general are 
BVIII and BIX. These models possess the highest dimension of space of structural con- 
stant. 

So, the integrable subclass of the Bianchi models forms a set of zero measure in the 
space of all Bianchi class A models. This fact means that integrable cases are excep- 
tional whereas those non-integrable ones are typical, as well as the Bianchi models which 
izotropize at infinite time in the full class of Bianchi models |1^. The establishing of the 



above fact for the Einstein equations in general is the challenge for the authors' future 
investigations. 
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